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Abstract 

In this paper, we investigate the transmission completion time minimization problem in a two-user 
additive white Gaussian noise (AWGN) broadcast channel, where the transmitter is able to harvest 
energy from the nature, using a rechargeable battery. The harvested energy is modeled to arrive at the 
transmitter randomly during the course of transmissions. The transmitter has a fixed number of packets 
to be delivered to each receiver. Our goal is to minimize the time by which all of the packets for 
both users are delivered to their respective destinations. To this end, we optimize the transmit powers 
and transmission rates intended for both users. We first analyze the structural properties of the optimal 
transmission policy. We prove that the optimal total transmit power has the same structure as the optimal 
single-user transmit power (TJ, J2). We also prove that there exists a cut-off power level for the stronger 
user. If the optimal total transmit power is lower than this cut-off level, all transmit power is allocated 
to the stronger user, and when the optimal total transmit power is larger than this cut-off level, all 
transmit power above this level is allocated to the weaker user. Based on these structural properties of 
the optimal policy, we propose an algorithm that yields the globally optimal off-line scheduling policy. 
Our algorithm is based on the idea of reducing the two-user broadcast channel problem into a single-user 
problem as much as possible. 

Index Terms 

Energy harvesting, rechargeable wireless networks, broadcast channels, transmission completion 
time minimization, throughput maximization. 

This work was supported by NSF Grants CCF 04-47613, CCF 05-14846, CNS 07-16311, CCF 07-29127, CNS 09-64632. 



1 



I. Introduction 

We consider a wireless communication network where users are able to harvest energy from 
the nature using rechargeable batteries. Such energy harvesting capabilities will make sustainable 
and environmentally friendly deployment of wireless communication networks possible. While 
energy-efficient scheduling policies have been well-investigated in traditional battery powered 
(un-rechargeable) systems S-jHl, energy-efficient scheduling in energy harvesting networks with 
nodes that have rechargeable batteries has only recently been considered (T), (2J. References 
EL ED consider a single-user communication system with an energy harvesting transmitter, and 
develop a packet scheduling scheme that minimizes the time by which all of the packets are 
delivered to the receiver. 

In this paper, we consider a multi-user extension of the work in JTJ, |0. In particular, we 
consider a wireless broadcast channel with an energy harvesting transmitter. As shown in Fig. \T\ 
we consider a broadcast channel with one transmitter and two receivers, where the transmitter 
node has three queues. The data queues store the data arrivals intended for the individual 
receivers, while the energy queue stores the energy harvested from the environment. Our objective 
is to adaptively change the transmission rates that go to both users according to the instantaneous 
data and energy queue sizes, such that the total transmission completion time is minimized. 

In this paper, we focus on finding the optimum off-line schedule, by assuming that the energy 
arrival profile at the transmitter is known ahead of time in an off-line manner, i.e., the energy 
harvesting times and the corresponding harvested energy amounts are known at time t = 0. We 
assume that there are a total of B 1 bits that need to be delivered to receiver 1, and B 2 bits 
that need to be delivered to receiver 2, available at the transmitter at time t = 0. As shown in 
Fig. [21 energy arrives (is harvested) at points in time marked with o; in particular, E k denotes 
the amount of energy harvested at time s k . Our goal is to develop a method of transmission to 
minimize the time, T, by which all of the data packets are delivered to their respective receivers. 

The optimal packet scheduling problem in a single-user energy harvesting communication 
system is investigated in JTJ, [|2]|. In JTJ, ||2), we prove that the optimal scheduling policy has a 
"majorization" structure, in that, the transmit power is kept constant between energy harvests, the 
sequence of transmit powers increases monotonically, and only changes at some of the energy 
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harvesting instances; when the transmit power changes, the energy constraint is tight, i.e., at the 
times when the transmit power changes, the total consumed energy equals the total harvested 
energy. In HI, 0, we develop an algorithm to obtain the optimal off-line scheduling policy 
based on these properties. Reference [9] extends fl}, to the case where rechargeable batteries 
have finite sizes. We extend [QQ|, m IflOl to a fading channel. 

References BU, IfTOll investigate two related problems. The first problem is to maximize the 
throughput (number of bits transmitted) with a given deadline constraint, and the second problem 
is to minimize the transmission completion time with a given number of bits to transmit. These 
two problems are "dual" to each other in the sense that, with a given energy arrival profile, if 
the maximum number of bits that can be sent by a deadline is B* in the first problem, then 
the minimum time to transmit B* bits in the second problem must be the deadline in the first 
problem, and the optimal transmission policies for these two problems must be identical. In 
this paper, we will follow this "dual problems" approach. We will first attack and solve the 
first problem to determine the structural properties of the optimal solution. We will then utilize 
these structural properties to develop an iterative algorithm for the second problem. Our iterative 
approach has the goal of reducing the two-user broadcast problem into a single-user problem as 
much as possible, and utilizing the single-user solution in UJ, 0. The second problem is also 
considered in the independent work ifTTl which uses convex optimization techniques to reduce 
the problem into local sub-problems that consider only two energy arrival epochs at a time. 

We first analyze the structural properties of the optimal policy for the first problem where 
our goal is to maximize the number of bits delivered to both users under a given deadline 
constraint. To that end, we first determine the maximum departure region with a given deadline 
constraint T. The maximum departure region is defined as the set of all (Bi,B 2 ) that can be 
transmitted to users reliably with a given deadline T. In order to do that, we consider the problem 
of maximizing [iiBi + H2B2 under the energy causality constraints for the transmitter, for all 
A*i, A*2 > 0. Varying /ii, /i 2 traces the boundary of the maximum departure region. We prove that 
the optimal total transmit power policy is independent of the values of /ii, /i 2 , and it has the 
same "majorization" structure as the single-user non-fading solution. As for the way of splitting 
the total transmit power between the two users, we prove that there exists a cut-off power level 
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for the stronger user, i.e., only the power above this cut-off power level is allocated to the weaker 
user. 

We then consider the second problem, where our goal is to minimize the time, T, by which a 
given (Bi, B 2 ) number of bits are delivered to their intended receivers. As discussed, since the 
second problem is "dual" to the first problem, the optimal transmission policy in this problem 
has the same structural properties as in the first problem. Therefore, in the second problem as 
well, there exists a cut-off power level. The problem then becomes that of finding an optimal 
cut-off power such that the transmission times for both users become identical and minimized. 
With these optimal structural properties, we develop an iterative algorithm that finds the optimal 
schedule efficiently. In particular, we first use the fact that the optimum total transmit power has 
the same structural properties as the single-user problem, to obtain the first optimal total power, 
Pi, i.e., the optimal total power in the first epoch. Then, given the fact that there exists a cut-off 
power level, P c , for the stronger user, the optimal transmit strategy depends on whether P 1 is 
smaller or larger than P c , which, at this point, is unknown. Therefore, we have two cases to 
consider. If P c is smaller than P\, then the stronger user will always have a constant, P c , portion 
of the total transmit power. This reduces the problem to a single-user problem for the second 
user, together with a fixed-point equation in a single variable (P c ) to be solved to ensure that 
the transmissions to both users end at the same time. On the other hand, if P c is larger than Pi, 
this means that all of Pi must be spent to transmit to the first (stronger) user. In this case, the 
number of bits delivered to the first user in this time duration can be subtracted from the total 
number of bits to be delivered to the first user, and the problem can be started anew with the 
updated number of bits (Bi, B' 2 ) after the first epoch. Therefore, in both cases, the broadcast 
channel problem is essentially reduced to single-user problems, and the approach in flTJ, [[2]] is 
utilized recursively to solve the overall problem. 

II. System Model and Problem Formulation 

The system model is as shown in Figs. \T\ and [2] The transmitter has an energy queue and two 
data queues (Fig. [[])• The physical layer is modeled as an AWGN broadcast channel, where the 
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received signals at the first and second receivers are 



Y 1 = X + Z 1 (1) 
Y 2 = X + Z 2 (2) 

where X is the transmit signal, and Z\ is a Gaussian noise with zero-mean and unit-variance, and 
Z 2 is a Gaussian noise with zero-mean and variance a 2 , where a 2 > 1. Therefore, the second user 
is the degraded (weaker) user in our broadcast channel. Assuming that the transmitter transmits 
with power P, the capacity region for this two-user AWGN broadcast channel is lfT2l 

n<ilog 2 (l + aP) (3) 

where a is the fraction of the total power spent for the message transmitted to the first user. 
Let us denote f(p) — |log 2 (1 +p) for future use. Then, the capacity region is n < f(aP), 
r 2 < / y ^p+S ) • Th^ s ca P ac ity region is shown in Fig. [3j 
Working on the boundary of the capacity region, we have 

P = 2 2 < ri+T *> + (a 2 - l)2 2r2 - a 2 (5) 
= 9(r u r 2 ) (6) 

As shown in Fig. \T\ the transmitter has B 1 bits to transmit to the first user, and B 2 bits to 
transmit to the second user. Energy is harvested at times Sk with amounts E^. Our goal is to 
select a transmission policy that minimizes the time, T, by which all of the bits are delivered to 
their intended receivers. The transmitter adapts its transmit power and the portions of the total 
transmit power used to transmit signals to the two users according to the available energy level 
and the remaining number of bits. The energy consumed must satisfy the causality constraints, 
i.e., at any given time t, the total amount of energy consumed up to time t must be less than or 
equal to the total amount of energy harvested up to time t. 

Before we proceed to give a formal definition of the optimization problem and propose the 
solution, we start with the "dual" problem of this transmission completion time minimization 
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problem, i.e., instead of trying to find the minimal T, we aim to identify the maximum number of 
bits the transmitter can deliver to both users by any fixed time T. As we will observe in the next 
section, solving the "dual" problem enables us to identify the optimal structural properties for 
both problems, and these properties eventually help us reduce the original problem into simple 
scenarios, which can be solved efficiently. 

III. Characterizing V{T): Largest B 2 ) Region for a Given T 

In this section, our goal is to characterize the maximum departure region for a given deadline 
T. We define it as follows. 

Definition 1 For any fixed transmission duration T, the maximum departure region, denoted as 
T>{T), is the union of (Bi, B 2 ) under any feasible rate allocation policy over the duration [0, T), 
i.e., V(T) = {j ri (t). r2 ( t ){Bi,B 2 ){ri{t),r2{t)), subject to the energy constraint f g(r 1 ,r 2 ){r)dr < 

Ei:s t <t E i>f 0r 0<t<T. 

We call any policy which achieves the boundary of V(T) to be optimal. In the single-user 
scenario in 0]|, we first examined the structural properties of the optimal policy. Based on these 
properties, we developed an algorithm to find the optimal scheduling policy. In this broadcast 
scenario also, we will first analyze the structural properties of the optimal policy, and then obtain 
the optimal solution based on these structural properties. The following lemma which was proved 
for a single-user problem in flTJ, [J2J was also proved for the broadcast problem in iflTTl . 

Lemma 1 Under the optimal policy, the transmission rate remains constant between energy 
harvests, i.e., the rate only potentially changes at an energy harvesting epoch. 

Proof: We prove this using the strict convexity of g(ri,r 2 ). If the transmission rate for any user 
changes between two energy harvesting epochs, then, we can always equalize the transmission 
rate over that duration without contradicting with the energy constraints. Based on the convexity 
of g(r 1 ,r 2 ), after equalization of rates, the energy consumed over that duration decreases, and 
the saved energy can be allocated to both users to increase the departures. Therefore, changing 
rates between energy harvests is sub-optimal. ■ 
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Therefore, in the following, we only consider policies where the rates are constant between 
any two consecutive energy arrivals. We denote the rates that go to both users as (r ln ,, r 2n ) over 
the duration [s n -i,s n ). With this property, an illustration of the maximum departure region is 
shown in Fig. HI 

Lemma 2 T>(T) is a convex region. 

Proof: Proving the convexity of V{T) is equivalent to proving that, given any two achievable 
points (£>i,£> 2 ) and (B[, B 2 ) in V(T), any point on the line between these two points is also 
achievable, i.e., in V{T). Assume that (fix, B 2 ) and (B[, B' 2 ) can be achieved with rate allocation 
policies (ri,r 2 ) and (r'^r^), respectively. Consider the policy (Ari + Ar / 1 ,Ar 2 + Ar 2 ), where 
A = 1 — A. Then, the energy consumed up to s n is 

n n n 

^g(\ru + \r' u , Xr 2i + \r' 2i )k < A ^ g(r u , r 2i )k + A ^g{r' u , r' 2i )k (7) 

i=l i=l i=l 

re— 1 re— 1 

<\Y,Ei+\J2 E * (8) 

i=0 i=0 
n-1 

Therefore, the energy causality constraint is satisfied for any A E [0, 1], and the new policy is 
energy-feasible. Any point on the line between (B±, B 2 ) and (B[, B' 2 ) can be achieved. When 
A 7^ 0, 1, the inequality in © is strict. Therefore, we save some amount of energy under the 
new policy, which can be used to increase the throughput for both users. This implies that V{T) 
is strictly convex. ■ 

In order to simplify the notation, in this section, for any given T, we assume that there are 
N — 1 energy arrival epochs (excluding t = 0) over (0,T). We denote the last energy arrival 
epoch before T as s N _i, and s N = T, with l N = T — s^-i, as shown in Fig. [5J 
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Since V{T) is a strictly convex region, its boundary can be characterized by solving the 
following optimization problem for all /ii,/i 2 > 0, 

TV N 

max /i! y~] r ln l n + fj, 2 Y] r 2n l n 

n=l n=l 

i j'-i 
s.t. ^g{r ln ,r 2n )l n <^E n , Vj:0<j<N (10) 

n=l n=0 

where l n is the length of the duration between two consecutive energy arrival instances s n and 
s n _i, i.e., l n — s n — s n -i, and ri and r 2 denote the rate sequences r\ n and r 2n for users 1 and 2, 
respectively. The problem in (flOl) is a convex optimization problem with a convex cost function 
and a convex constraint set, therefore, the unique global solution should satisfy the extended 
KKT conditions. 
The Lagrangian is 

N N 



n=l n=l 
N / j j-l \ N N 



£(ri, r 2 , A, 7) =fl 1 ^ r lJn + /X 2 ^ 
n=l n=l 
/V / j 

— Aj- I ^2g(r ln ,r 2n )l n ~^2 E n\ + J^7ln r ln + ^ 72n^2n (H) 

j=l \n=l n=0 / n=l n=l 

Taking the derivatives with respect to r ln and r 2n , and setting them to zero, we have 

Mi + 7m - ( X; A,) 2 2 ^"+^) = 0, n = l,...,N (12) 

\j"=n / 

/i 2 + 72n - ( X i) (2 2(ri " +r2 " } + (a 2 - l)2 2r -) =0, n = 1, . . . , N (13) 

\j=n / 

where 7i n = if r± n > 0, and 7 2n = if r 2n > 0. Based on these KKT optimality conditions, 
we first prove an important property of the optimal policy. 

Lemma 3 The optimal total transmit power of the transmitter is independent of the values of 
f/,1, /i 2 , and it is the same as the single-user optimal transmit power. Specifically, 

> (14) 



P n = J ~ t "" 1 (15) 



i.e., at t = si n , P n switches to P n +\- 

Proof: Based on the expression of g(r ln ,r 2n ) in © and the KKT conditions in (fT2~l)-(fT3l). we 
have 

g{r ln , r 2n ) = - a (16) 

> 2 2 ( ri ™+ r2 ™) — i (17) 

where (fTTT ) becomes an equality when r 2n = 0. Therefore, when r 2n > 0, ([T6l)-([T9l imply 

9(rm, r 2n ) = -g— - a 2 > - 1 (20) 

When r 2n = 0, we must have r ln > 0. Otherwise, if r ln = 0, we can always let the weaker user 
transmit with some power over this duration without contradicting with any energy constraints. 
Since there is no interference from the stronger user, the departure from the weaker user can 
be improved, thus it contradicts with the optimality of the policy. Therefore, when r 2 ,„ = 0, 
7i„ = 0, and (fT6l)-(fl~9l imply 

g(r ln , r 2n ) = ^ - 1 > ^ 2 - a 2 (21) 

Therefore, we can express g(r ln ,r 2n ) in the following way: 

g(ri n , r 2n ) = max <^ jf 1 - 1, ^ - a 2 I (22) 

Plotting these two curves in Fig. [61 we note that the optimal transmit power, P n = g(r ln ,r 2n ), 
is always the curve on the top. If rf' 2 a 2 > tf 1 1 for some n, then, we have 

vj v : > vj v : > ° ~ !. Vri > n ( 23 ) 
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where the first inequality follows from the KKT condition that Xj > for j = 1,2, ... N. 
Therefore, we conclude that there exists an integer n, < n < N, such that, when n < n, 
r 2n = 0; and when n > n, r 2n > 0. 

Furthermore, (f20l)-(f2TT) imply that, the energy constraint at t — s n must be tight. Otherwise, 
X n = 0, and (12TI) implies 

g{rm, r 2 n) = N ^ ~ 1 > ^ - a 1 = g(r hn+1 , r 2 ,n+i) (24) 



which contradicts with (1201) . Therefore, in the following, when we consider the energy constraints, 
we only need to consider two segments [0, s n ) and [s n+ i, sjv) separately. 

When n < n, based on (l20l) . if A n = 0, we have g(r ln ,r 2n ) = g(r ljn+ i, r 2>n+ i). Starting 
from n — 1, g(ri n ,r 2n ) remains a constant until an energy constraint becomes tight. Therefore, 
between any two consecutive epochs, when the energy constraints are tight, the power level 
remains constant. Similar arguments hold when n > n. Thus, the corresponding power level is 

Ein — 1 771 
i-i 1 
F n = i- 1 - 1 (25) 

S in ~ S in-1 

where s in _ 1 and s in are two consecutive epochs with tight energy constraint. 

Finally, we need to determine the epochs when the energy constraint becomes tight. Another 
observation is that g(r m ,r 2n ) must monotonically increase in n, as shown in Fig. [6] This is 
because both of these two curves monotonically increase, and the maximum value of these two 
curves should monotonically increase also. Therefore, based on the monotonicity of the transmit 
power, we conclude that 

d = i r, 

i„-i<i<N I S; — Si 



i n = arg min I — J %n 1 \ (26) 



'ln-1 



This completes the proof. ■ 

Since the power can be obtained directly irrespective of the values of /ii, fi 2 , the optimization 
problem in (flOl) is separable over each duration [s n _i, s n ). Specifically, for < n < N, the local 
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optimization becomes 



max H\Tx n + /i 2 r 2n 

rin,T2n 

s.t. g(r ln ,r 2n ) < P n (27) 

We relax the power constraint to be an inequality to make the constraint set convex. Thus, this 
becomes a convex optimization problem. This does not affect the solution since the objective 
function is always maximized on the boundary of its constraint set, i.e., the capacity region 
defined by the transmit power P n . 

When — < o n t\ » me solution to (1271 can be expressed as 

r ln = Uog 2 {l + P n ) (28) 
r 2 „ = (29) 

In this scenario, all of the power P n is allocated to the first user. 
When < & < a 2 , we have 

nn=hoJ^^) (30) 



2 V A*2 - A*i 

1, /(/i 2 -/ii)(P n + (7 2> 

2 V /W 2 - 1 



r 2 „ = - log 2 ( — 2 ) (31) 



In this scenario, a constant amount of power, — 1, is allocated to the first user, and the 

remaining power is allocated to the second user. 

When ^ > cr 2 , we have 
mi 

Tin = (32) 
Pn 

2" oz V" ' (J 2 



r 2n = -log 2 (l + ^) (33) 



In this scenario, all of the P n is allocated to the second user. 
Let us define a constant power level as 



1^2 ~ fJ>l 
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Based on the solution of the local optimization problem in (|27T) . we establish another important 
property of the optimal policy as follows. 

Lemma 4 For fixed ji 2 , under the optimal power policy, there exists a constant cut-off power 
level, P c , for the first user. If the total power level is below this cut-off power level, then, all the 
power is allocated to the first user; if the total power level is higher than this level, then, all 
the power above this cut-off level is allocated to the second user. 

In the proof of Lemma |3j we note that the optimal power P n monotonically increases in n. 
Combining Lemma |3] and Lemma HI we illustrate the structure of the optimal policy in Fig. [7J 
Moreover, the optimal way of splitting the power in each epoch is such that both users' shares of 
the power monotonically increase in time. In particular, the second user's share is monotonically 
increasing in time. Hence, the path followed in the (J3i, B 2 ) plane is such that it changes direction 
to get closer to the second user's departure axis as shown in Fig. @] The dotted trajectory cannot 
be optimal, since the path does not get closer to the second user's departure axis in the middle 
(second) power epoch. 

Corollary 1 Under the optimal policy, the transmission rate for the first user, {r ln }^ =1 , is 
either a constant sequence (zero or a positive constant), or an increasing sequence. Moreover, 
before r\ n achieves its final constant value, r 2n = 0; and when n„ becomes a constant, r 2n 
monotonically increases in n. 

Based on Lemma |3j we observe that for fixed T, fix and fi 2 , the optimal total power allocation 
is unique, i.e., does not depend on {i x and [i 2 . However, the way the total power is split between 
the two users depends on fi 2 . In fact, the cut-off power level P c varies depending on the 
value of fi 2 /fi±. Therefore, for different values of fi 2 / ' Hi, the optimal policy achieves different 
boundary points on the maximum departure region, and varying the value of fi 2 ] 7*i traces the 
boundary of this region. 

In this section, we characterized the maximum departure region for any given time T. We 
proved that the optimal total transmit power is the same as in the single-user case, and there 
exists a cut-off power for splitting the total transmit power to both users. In the next section, we 
will use these structural properties to solve the transmission completion minimization problem. 
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IV. Minimizing the Transmission Completion Time T for a Given (B 1 , B 2 ) 

In this section, our goal is to minimize the transmission completion time of both users for a 
given (Bx, B 2 ). The optimization problem can be formulated as 



mm 

ri,r 2 



T 



S.t. 



3 J'-l 



Vj : < j < N(T) 



n=l n=l 




^2 r lnL =B X , ^ T ^il 



n 



B 2 



(35) 



n=l n=l 



where N(T) — 1 is the number of energy arrival epochs (excluding t = 0) over (0,T), and 
In(t) = T — sn(t)-i- Since N(T) depends on T, the optimization problem in (1351) is not a 
convex optimization problem in general. Therefore, we cannot solve it using standard convex 
optimization tools. 

We first note that this is exactly the "dual" problem of maximizing the departure region for 
fixed T. They are "dual" in the sense that, if the minimal transmission completion time for 
(Bx,B 2 ) is T, then (B t , B 2 ) must lie on the boundary of V(T), and the transmission policy 
should be exactly the same for some (jji,fx 2 ). This is based on the fact the V(T) C V{T') for 
any T < T'. Assume (B 1 ,B 2 ) does not lie on the boundary of V(T). Then, either (B 1: B 2 ) 
cannot be achieved by T or (Bi, B 2 ) is strictly inside V{T) and hence (Bi, B 2 ) can be achieved 
by T' < T. Therefore, if (Bx, B 2 ) does not lie on the boundary of V(T), then T cannot be the 
minimum transmission completion time. 

We have the following lemma. 

Lemma 5 When Bx,B 2 ^ 0, under the optimal policy, the transmissions to both users must be 
finished at the same time. 

Proof: This lemma can be proved based on Corollary \T\ If the transmission completion time for 
both users is not the same, then over the last duration, we transmit only to one of the users, 
while the transmission rate to the other user is zero. This contradicts with the monotonicity of 
the transmission rates for both users. Therefore, under the optimal policy, the transmitter must 
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finish transmitting to both users at the same time. ■ 

This lemma is proved in [fTTTl also, by using a different approach. The authors prove it in IfTTI 
mainly based on the convexity of the capacity region of the broadcast channel. 

For fixed (B\,B 2 ), let us denote the transmission completion time for the first and second 
user, by T\ and T 2 , respectively. We note that T\ and T 2 depend on the selection of the cut-off 
power level, P c . In particular, T\ is monotonically decreasing in P c , and T 2 is monotonically 
increasing in P c . Based on Lemma [5l the problem of optimal selection of P c , can be viewed 
as solving a fixed point equation. In particular, P c must be chosen such that, the resulting 7\ 
equals T 2 . Therefore, we propose the following algorithm to solve the transmission completion 
time, T, minimization problem. Our basic idea is to try to identify the cut-off power level P c in 
an efficient way. 

Since the power allocation is similar to the single-user case (c.f. Lemma [3]), our approach 
to find T will be similar to the method in [0Q, fl2]. First, we aim to identify P x , the first total 
transmit power starting from t = in the system. This is exactly the same as identification of 
Pi in the corresponding single-user problem. For this, as in (T), 0, we treat the energy arrivals 
as if they have arrived at time t — 0, and obtain a lower bound for the transmission completion 
time as in Q], ||2]. In order to do that, first, we compute the amount of energy required to finish 
(£>i, B 2 ) by s\. This is equal to g s\, denoted as A\. Then, we compare A\ with E . If 

E is greater than A x , this implies that the transmitter can finish the transmission before s\ with 
E , and future energy arrivals are not needed. In this case, the minimum transmission completion 
time is the solution of the following equation 

g (f , f ) T = E (36) 

If Ai is greater than E , this implies that the final transmission completion time is greater than 
si, and some of the future energy arrivals must be utilized to complete the transmission. We 
calculate the amount of energy required to finish (B x , B 2 ) by s 2 , s 3 , . . . , and denote them as A 2 , 
A 3 , . . ., and compare these with E + E 1 , Y^j=o Z^j=o Ej, ■ ■ ■ , until the first A+ that becomes 
smaller than Y^=o ^ e denote the corresponding time index as i\. Then, we assume that we 
can use X^Lo* t0 transmit (Bi,B 2 ) at a constant rate. And, the corresponding transmission 
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completion time is the solution of the following equation 

,(**)r = X> (37) 

We denote the solution to this equation as T, and the corresponding power as Pi. From 
our analysis, we know that the solution to this equation is the minimum possible transmission 
completion time we can achieve. Then, we check whether this constant power P x is feasible, 
when the actual energy arrival times are imposed. If it is feasible, it gives us the minimal 
transmission completion time; otherwise, we get P 1 by selecting the minimal slope according 
to (fT5l) . That is to say, we draw all of the lines from t = to the corner points of the energy 
arrival instances before T, and choose the line with the smallest slope. We denote by the 
corresponding duration associated with P±. This is shown in Fig. [8] 

Once Pi is selected, we know that it is the optimal total transmit power in our broadcast 
channel problem. Next, we need to divide this total power between the signals transmitted to 
the two users. Based on Lemma @] and Corollary [T} if the cut-off power level P c is higher than 
Pi, then, the transmitter spends all Pi for the stronger user; otherwise, the first user finishes its 
transmission with a constant power P c . 

We will first determine whether P c lies in [0, Pi] or it is higher than Pi. Assume P c = Pi. 
Therefore, the transmission completion time for the first (stronger) user is 

r, = ^ (38, 

Once P c is fixed, we can obtain the minimum transmission completion time for the second user, 
T 2 , by subtracting the energy consumed by the first user, and treating Pi as an interference for 
the second user. This reduces the problem to the single-user problem for the second user with 
fading, where the fading level is Pi + <r 2 over [0, Ti), and a 2 afterwards. The single-user problem 
with fading is studied in IfTOll . Since obtaining the minimal transmission completion time is not 
as straightforward for the fading channel, a more approachable way is to calculate the maximum 
number of bits departed from the second user by Ti, denoted as D 2 (Ti, P c ). In order to do that, 
we first identify the optimal power allocation policy with fixed deadline Ti. This can be done 
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according to Lemma [3] Assume that the optimal power allocation gives us Pi, P 2 , ■ ■ ■ , fjv(Ti)- 
Then, we allocate Pi to the first user over the whole duration, and allocate the remaining power 
to the second user. Based on (|4]), we calculate the transmission rate for the second user over 
each duration, and obtain D 2 (Ti,P c ) according to 

1 / P -P\ 
D 2 (T h P c ) = g - log \l + (s in - Sin _J (39) 

We observe that, given P c , D 2 (Ti,P c ) is a monotonically increasing function of Ti. Moreover, 
given Ti, D 2 (Ti,P c ) is a monotonically decreasing function of P c . 

If D 2 (Ti, P c ) is smaller than B 2 , it implies that Ti < T 2 , and we need to decrease the rate for 
the first user to make Ti and T 2 equal. Based on Lemma 0] this also implies that the transmission 
power for the first user is a constant P c < P\. In particular, P c is the unique solution of 

Note that D 2 (j^|ry> Pcj ls a continuous, strictly monotonically decreasing function of P c , hence 
the solution for P c in (|4Q|) is unique. Since Ti is a decreasing function of P c and _D 2 (/(p - )) Pcj 
is a decreasing function of P c , we can use the bisection method to solve (l40l) . In this case, the 
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minimum transmission completion time is T - 

If D 2 (Ti, P c ) is larger than B 2 , that implies T 2 < T 1? and we need to increase the power 
allocated for the first user to make Ti and T 2 equal, i.e., P c > Pi. Therefore, from Lemma H 
over the duration [0, SfJ, the optimal policy is to allocate the entire Pi to the first user only. We 
allocate Pi to the first user, calculate the number of bits departed for the first user, and remove 
them from B>i. This simply reduces the problem to that of transmitting (B[, B> 2 ) bits starting at 
time t = Sj 15 where B[ = B x — /(P^s^. The process is illustrated in Fig. [91 Then, the minimum 
transmission completion time is 

r = s « + 1(E) <41) 

where K is the number of recursions needed to get P c . 

In both scenarios, we reduce the problem into a simple form, and obtain the final optimal 



16 



policy. Before we proceed to prove the optimality of the algorithm, we introduce the following 
lemma first, which is useful in the proof of the optimality of the algorithm. 

Lemma 6 f(E/T)T monotonically increases in T; f {jjz^§ffy^s^ T monotonically increases 
in T also. 

Proof: The monotonicity of both functions can be verified by taking derivatives, 

(f(E/T)Ty = f(E/T)- (2]n2 * T + E) (42) 

and 

U(EIT)T)» = JL - ^L_) < (43) 

where the inequality follows since E > 0. Therefore, f(E/T)T is a strictly concave function, and 
its first derivative monotonically decreases when T increases. Since when lim T _^ 00 (/(_E/T)T) / 
= 0, when T < oo, we have (f(E/T)T)' > 0, therefore, the monotonicity follows. 
Similarly, we have 

(' ( (i-S^ ) T ) =5 log2 W ' + E/T) - 5 ' og2 ( ° 2 + (1 " a)EIT) 

EE E (l-a)E 

+ 01 ^ ^ (44) 



21n2£ + (T 2 r 21n2(l - a)E + a 2 T 

and 



/f aE/T ) T Y = _&_( 1 ^^ <0 (45) 

J \(l-aE/t) + a 2 J J 2T\n2\{a 2 T/{l-a) + E) 2 (a 2 T + E) 2 ) V ; 

Again, the concavity implies that the first derivative is positive when T < oo, and the mono- 
tonicity follows. ■ 

Theorem 1 The algorithm is feasible and optimal. 

Proof: We first prove the optimality. In order to prove that the algorithm is optimal, we need 
to prove that Pi is optimal. Once we prove the optimality of P x , the optimality of P 2 , P 3 , 
. . . follows. Since the solution obtained using our algorithm always has the optimal structure 
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described in Lemma 0] the optimality of the power allocation also implies the optimality of 
the rate selection, thus, the optimality of the algorithm follows. Therefore, in the following, we 
prove that Pi is optimal. 

First, we note that Pi is the minimal slope up to T. We need to prove that Pi is also the 
minimal slope up to the final transmission completion time, T. Let us define T' as follows 

Pi 

Assume that with Pi, we allocate aPi to the first user, and finish (Pi, P 2 ) using constant rates. 
Then, we allocate aPi to the first user, and the rest to the second user. Based on Lemma [6l we 
have 

/(aPi)T' > /(aA)f = B, (47) 

/( ,, °* i )r>f( n 2 )t = B 2 (48) 

\{l-a)Pi + a 2 J y(\ - a ) Pi + a 2 ) 

Therefore, T' is an upper bound for the optimal transmission completion time. Since Pi is the 
minimal slope up to T', we conclude that Pi is optimal throughout the transmission. Following 
similar arguments, we can prove the optimality of the rest of the power allocations. This 
completes the proof of optimality. 

In order to prove that the allocation is feasible, we need to show that the power allocation for 
the first user is always feasible in each step. Therefore, in the following, we first prove that Pi 
is feasible when we assume that P c — P\. The feasibility of Pi also implies the feasibility of 
the rest of the power allocation. With the assumption that P c = P\, the final transmission time 
for the first user is 

Ti = < Bl (49) 

1 /(Pi) " f{aP x ) 1 V 

Based on (HTJl and (|48l> . we know that T x < V . Since Pi is feasible up to V, therefore, Pi 
is feasible when we assume that P c = Pi. The feasibility of the rest of the power allocations 
follows in a similar way. This completes the feasibility part of the proof. ■ 
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V. Numerical Examples 

We consider a band-limited AWGN broadcast channel, with bandwidth W = 1 MHz and 
the noise power spectral density iV = 10~ 19 W/Hz. We assume that the path loss between the 
transmitter and the first receiver is about 100 dB, and the path loss between the transmitter and 
the second user is about 105 dB. Then, we have 

ri = w l0fe t 1 + w) = log2 + 5p) Mbps <50) 

Therefore, 

^( ri , r 2 ) = 10~ 3 2 ri+r2 + (10~ 2 ' 5 - 10" 3 )2 ra - 10~ 2 ' 5 W (52) 

For the energy harvesting process, we assume that at times t = [0,2,5,6,8,9,11] s, we have 
energy harvested with amounts E = [10, 5, 10, 5, 10, 10, 10] mJ. We find the maximum departure 
region V{T) for T = 6, 8, 9, 10 s, and plot them in Fig. [101 We observe that the maximum 
departure region is convex for each value of T, and as T increases, the maximum departure 
region monotonically expands. 

Then, we aim to minimize the transmission completion time with (Bi,B 2 ) = (15,6) Mbits. 
Following our algorithm, we obtain the optimal transmission policy, which is shown in Fig. [TT] 
We note that the powers change only potentially at instances when energy arrives (Lemma [D; 
power sequence is monotonically increasing and "majorized" over the whole transmission dura- 
tion (Lemma [3]). We also note that, for this case, the first user transmits at a constant rate, and 
the rate for the second user monotonically increases. The transmitter finishes its transmissions 
to both users by time T = 9.66 s, and the last energy harvest at time t = 11 s is not used. 

Next, we consider the example when (Bi,B 2 ) = (20,2) Mbits, we have the optimal trans- 
mission policy, as shown in Fig. [121 In this example, the cut-off power is greater than Pi, and 
therefore, Pi is allocated to the first user only over [0, 5) s, and after t — 5 s, the first user 
keeps transmitting at a constant rate until all bits are transmitted. In this case, the transmission 
rates for both users monotonically increase. The transmitter finishes its transmissions by time 
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T = 9.25 s, and the last energy harvest is not used. 

VI. Conclusions 

We investigated the transmission completion time minimization problem in an energy harvest- 
ing broadcast channel. We first analyzed the structural properties of the optimal transmission 
policy, and proved that the optimal total transmit power has the same structure as in the single- 
user channel. We also proved that there exists a cut-off power for the stronger user. If the optimal 
total transmit power is lower than this cut-off level, all power is allocated to the stronger user, 
and when the optimal total transmit power is greater than this cut-off level, all power above this 
level is allocated to the weaker user. Based on these structural properties of the optimal policy, 
we developed an iterative algorithm to obtain the globally optimal off-line transmission policy. 
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Fig. 1. An energy harvesting two-user broadcast channel. 



O RX 1 



O RX 2 





b 




E, 

-J - 


E K 

i 




'0 


•Si 


S2 


Sk 


T 


t 



{B X ,B 2 ) 

Fig. 2. System model. (Bi, B2) bits to be transmitted to users are available at the transmitter at the beginning. Energies arrive 
(are harvested) at points denoted by o. T denotes the transmission completion time by which all of the bits are delivered to 
their respective destinations. 




Ci Ri 



Fig. 3. The capacity region of the two-user AWGN broadcast channel. 
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Fig. 4. The maximum departure region and trajectories to reach the boundary. Dotted trajectory is not possible. 
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Fig. 5. Rates (rim ran) an d corresponding durations l n with a given deadline T. 
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Fig. 6. The value of the optimal transmit power is always equal to the curve on top. 
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Fig. 7. Optimally splitting the total power between the signals that go to the two users. 
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Fig. 8. Determining the optimal total power level of the first epoch. 
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Fig. 9. Search for the cut-off power level P c iteratively. 
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Fig. 10. The maximum departure region of the broadcast channel for various T. 
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Fig. 11. Cut-off power P c = 1.933 mW. Optimal transmit rates are ri = 1.552 Mbps, r 2 = [0.274, 0.680, 1.369, 1.834] 
Mbps, with durations 1 = [5,3, 1,0.66] s. 
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Fig. 12. Cut-off power P c = 4.107 mW. Optimal transmit rates n = [2,2.353,2.353,2.353] Mbps and r 2 = 
[0, 0.167, 0.856, 2.570] Mbps, with durations 1 = [5, 3, 1, 0.25] s. 
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